Quantum Phase Transition in Skyrmion Lattice 
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We investigate the ground state of 2D electron gas in Quantum Hall regime at the filling factor 
' slightly deviating from unity, that can be viewed as a sparse lattice of skyrmions. We have found 

y—( | that in the low density limit skyrmions are bound in pairs, those forming the actual lattice. We 

, | . have shown that at further density increase the lattice undergoes a quantum phase transition, an 

analogue of superconducting phase transition in Josephson junction arrays. 

^ ! PACS numbers: 73.40.Hm; 73.20.Dx; 12.39.Dc. 

\ Skyrmion [jjj has been a favorite model particle of field theorists for decades. Only recently it has been comprehended 
that the low- lying charged excitations of the Integer Quantum Hall ground state may be skyrmions ( see || for 
' j review). This interesting development has prompted intensive experimental studies that have proved these exotic 
i particles to be real indeed. A skyrmion can be viewed as a topologically non-trivial distortion of the spins of many 
ground state electrons. That makes its spin much bigger than unity, although the skyrmion bears a unitary charge. 
^ The size of the skyrmions is controlled by a parameter g, the ratio of Zeeman energy to exchange energy per particle. 
D ' The smaller the g is, the bigger is the skyrmion. 

If the filling factor slightly deviates from unity, the low-lying charged excitations must appear at the background 
of the Quantum Hall state to compensate for excess charge. Therefore, the deviation of the filling factor sets the 
d . skyrmion concentration and the many-skyrmion ground state can be easily made by changing either magnetic field or 
electron density of 2D gas. The two studies of many-skyrmion ground state at finite g have been recently reported. 
HQ Despite the different methods applied and the different results obtained, it has been assumed in both studies that 
the skyrmions form a plain regular lattice. Also, in both studies no quantum effects have been taken into account. 

Our results demonstrate that the formation and development of the multi-skyrmion state is more complicated and 
interesting than it has been assumed previously. We have shown that at low concentrations the skyrmions always 
appear in pairs. These pairs form a triangular Wigner lattice. With increasing concentration the distance between 
I ' pairs decreases, and the interaction between the pairs strengthens. This leads to a quantum transition in skyrmion 
;> lattice. Below the transition point, the z-component of spin of each pair is well-defined. We call the variable that 
0^ ' is canonically conjugated to z-component of spin the phase. Above the transition point, it is this phase that is 
, well-defined. 

' We begin with the energy functional that describes spin textures in the limit of small g (see e.g. 
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It contains three terms that describe respectively spin stiffness, Zeeman energy, and electrostatic energy. It depends 
on m(r), a unitary vector characterizing local spin density. We adopt the microscopic expression S for the spin 
stifness at filling factor 1: p s = (l/16-\/27r)(e 2 /eA), A being the magnetic length. In the Zeeman term, g stands for the 
electron g-factor, and hb is the Bohr magneton. In the electrostatic term, q(r) presents the deviation of the electron 
density from the uniform background density n/, = l/27rA 2 , q being related to the density of topological charge B: 



-^e M „(m[<9 M m x <9„m]). (2) 
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5-H ' The precise definition of g reads g = g^sB / ' E ex , E ex = e /e\. 

The stiffness term dominates the energy and must be minimized first of all. It has been shown by Belavin and 
Polyakov Q that the textures minimizing the stiffness term are highly degenerate. They parametrized the spin density 
as follows: 

w m x + im y 2W \W\ 2 -1 

It turns out that all the textures corresponding to any analytical function W{z = x + iy) with N poles possess the 
same energy N \fn /32E ex . Neither pole (skyrmion positions) nor pole residuals (skyrmion radii) are fixed. Those are 
determined by the interplay of the weaker interactions: Zeeman and electrostatic energy. 
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Let us start with a single skyrmion. m We take W(z) — a/z, a being the skyrmion radius to be evaluated. 
The Coulomb energy is given by E Q = (37r 2 /64)e 2 /ea. The integral that gives Zeeman energy appears to diverge 
logarithmically at long distances from the core, since x, y spin components decrease very slowly with increasing 
distance. A slightly improved calculation shows that the divergency is to be cut at r s ~ A<7 -1 / 2 , and Zeeman energy 
reads Ez = gfiBB(a 2 / X 2 )lu.(r s /a) provided a -C r s . The minimization of the total energy with respect to a yields 
a ~ A(7n(l/g)g) -1 / 3 , so that a <C r s indeed. 

Let us consider now two skyrmions with opposite phases, 

. . a a 

W(z) = — -, 4 

z + R z — R 

the minus sign before the second term accounting for that phase shift. We see that the x, y spin components are 
effectively quenched at \z\ 3> R, that lowers Zeeman energy. For well-separated skyrmions (a <C R) the logarithmical 
divergency is cut at distances of the order of R and Zeeman energy reads 2gfiBB(a 2 / X 2 ) ln(R/a). This results in an 
attractive force between skyrmions oc 1/R. At long distances that always prevails Coulomb repulsion, so that the 
skyrmions are bound to form pairs. To find how they are exactly bound, we have studied the question numerically. 
Fig. 1 shows our results for the two skyrmion energy as function of their separation R, at optimal a. It appears that 
the minimum is achieved at zero separation, corresponding to W(z) = b 2 /z 2 . Two skyrmions merge as it sketched 
in Fig. 1. The optimal value of \b\ is bo ~ A/5 1 / 3 . The skyrmion pair bears big spin, So — 0.78/g 2 / 3 . We checked 
numerically that skyrmion attraction saturates, that is, there is no further bounding of pairs. 

Therefore we have shown that the charged excitations with lowest energy are eventually skyrmion pairs bearing 
charge ±2e. We proceed now with analysing many-skyrmion ground state. |L^ ] Let us compare now Zeeman and 
Coulomb contributions to the interaction between pairs. The latter appears to prevail provided the separation between 
pairs exceeds their equilibrium radius bo. Therefore, the pairs behave very much like point-like charges and form a 
triangular Wigner lattice to minimize their Coulomb energy. Such an arrangement persist up to relatively high 
skyrmion densities n s k y /n e i ~ (60/ A) -2 ~ g 2 ^ where they begin to overlap. Most probably, the overlap would depair 
skyrmions and one of the single-skyrmion lattices proposed would be realized at higher densities. In the present paper, 
we restrict our attention to lower densities where the pairs certainly form the triangular lattice. 

To reveal the interesting physics that persists even in this limiting case, let us note that the two-skyrmion solution 
found is infinitely degenerate. All the spins can be rotated about z-axis by an arbitrary angle yielding the texture 
of the same energy. Or, alternatively, we can multiply W(z) by an arbitrary phase factor exp(i(p). We will call the 
variable parametrizing all the degenerate textures the phase. This degeneracy is not physical and arises from the fact 
that we treat the pair, which is a quantum particle, as a classical object. 

Let us consider quantization of an isolated pair. We note that the shift in the phase space, <f> — ► <fi + 9, is just 
a rotation in spin space about z-axis being represented by the quantummechanical operator exp(— i9S z ). Here S z 
stands for the total spin of the skyrmion. Expanding near 9 = 0, we find that S z — —id/d<f>. That leads us to an 
important conslusion: 

(5) 

the phase and the total spin of the pair are conjugate variables satisfying Heisenberg uncertainty relation. Owing to 
degeneracy with respect to phase, physical isolated skyrmion pairs do not have a certain phase but rather possess a 
certain integer z-component of spin. This picture has very much in common with Coulomb blockade in superconducting 
island [ pT|JlC| ], that does not have a certain superconducting phase but rather a certain integer charge. 

Since the total spin of the pair is a function of its radius b, as well as its energy, the spin quantization leads to the 
energy quantization of pair states. The spectrum can be obtained from classical dependences E(b) and S z (b). Near 
the minimum energy it reads 

E = E Q + E c (S z -So) 2 , E c ^-~^. (6) 

Here Eq , Sq are optimal classical energy and spin respectively. Let us note that Sq is a continious function of g, to be 
contrasted with integer S z . If we change g, Sq changes continuously, whereas the optimal S z jumps between integer 
values. In the point of jump, energies of two states with different S z match providing extra degeneracy. 

This is reminiscent of Coulomb blockade phenomenon in Josephson arrays. JlC| ] Josephson arrays consist of super- 
conducting islands, the number of discrete charges in each island being canonically conjugate to its superconducting 
phase. The actual state of the array is determined by an interplay of charging energy and Josephson coupling. If 
charging energy dominates, the array is in insulating phase. There is a gap for charged excitations: one has to pay 
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charging energy in order to add or extract a discrete charge to/from the array. If the Josephson coupling increases, 
the quantum phase transition to the superconducting phase takes place. The long-range phase ordering appears, the 
gap vanishes. 

We shall expect similar behaviour from skyrmion lattices, with discrete charge replaced by discrete S z . At very 
low skyrmion concentrations the pairs can be regarded as isolated ones each having a well-defined spin and strongly 
fluctuating phase. The x, y components of the local spin density are undefined and have zero average, as it is sketched 
in Fig. 2a. There is a "spin gap" in this " inslulating" phase: one has to pay energy of the order of E c to increase 
or decrease the spin of the system by unity. If we increase skyrmion concentration, the pairs come closer to each 
other. That increases coupling, that tries to fix the phases of neighboring pairs. Therefore we expect a quantum 
transition to "superconducting" phase. At the transition point, the "spin gap" vanishes to zero. Instead, long-range 
phase ordering appears. That is, x, y components of the local spin density become non-zero. 

To give a quantitative description, we evaluate first the phase-dependent interaction between pairs. It appears that 
this interaction measures up E c at interskyrmion distances much larger than screening radius r s , (r s /A) 2 = \JH j (4v2<?) 
this is why we shall evaluate it for interpair distances r in the region r» r s , b. We expand energy functional up to 
the terms quadratic in s XjV and minimize the resulting expression matching s x . v near the pair cores with asymptotics 
of W(z). That yields 
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E int = £jcos(0i - fa), Ej = gf i B B^exp{ )(-£y/irr,r/2). (7) 

r s r s A 

As a result, we have the following many-body Hamiltonian, which takes into account "charging" and Josephson 
energies: 

H = E E cH^- - So) 2 + 1 -±E J cosfa - fr), (8) 



'j 



where E c is given by Eq.(^), Ej > is given by Eq.(0). Indexes label sites of triangular lattice, the prime restricts 
sum over j to the nearest neighbours of a certain site. Note that in our case the "Josephson" coupling has an 
antiferromagnetic sign, to be contrasted with Josepson arrays. In the limit Ej ^> E c the phases become well-defined 
and we can neglect "charging" term. The Hamiltonian reduces to energy functional. In the ground state, the lattice 
subdivides into three sublattices. The phases of sublattices are rotated by 2tt/3 with respect to each other (Fig 
2b). 

We are interested in characterizing the transition. Since the exact solution of the Hamiltonian is beyond the reach, 
we find the transition point in Hartree approximation. We approximate the exact ground-state wave function by a 
product of site wave functions, 

N 

*(0i,...,0iO=n*<(&)> ^ 

i=l 

to be determined from the minimum energy principle. The site wave functions then satisfy the following mean-field 
Hamiltonian: 

d 1 ' 

H M f = E Ec ^~ i Jhh ~ S °^ 2 + 7 Ej ^2^ n 3 ex P (-*</>*) - n i n * 3 + c - c -\- 

i 1 ij 

where we introduce the order parameter 
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at each site. We assume that near transition point the order parameter has the same symmetry as for Ej 3> E c . The 
non-trivial solution for the order parameter first appear at 

E°j = ^-{l-4(lS + 1/2] -So) 2 }, (10) 

this equation determining the transition point. Making use of Eqs.(|^, |?], ^) we find that the transition occurs at 
skyrmion density n s k y ~ rj 2 In - (!/</)■ To give more detailed predictions, we plot the critical skyrmion density versus 
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g in realistic range of this parameter. The curve exhibits sharp drops each time the energies of two spin states match, 
since this favours spin fluctuations and, consequently, facilitates phase ordering. 

To observe the transition experimentally, one shall measure either spin gap or the local magnetization parallel to the 
plane, those quantities respectively vanish and appear above transition point. Both quantities are likely to be accessed 
in magnetic resonance measurement. The experiment shall be performed at temperature that is lower than the typical 
energy scale involved, E c . Our estimations for n e i — 10 16 m~ 2 , g = 0.025 give E c /kB — IK. The fine tuning of the 
filling factor in the range of 0.01 is required. Our model does not explicitly include disorder. However we expect that 
the quantum transition occurs even in disordered heterostructures, provided the skyrmions survive. The disorder may 
trap the skyrmions in random potential minima, thus preventing the formation of any regular lattice. Albeit it can 
not pin the skyrmion phase, that to be pinned only by interskyrmion interactions with increasing concentration. This 
is why we expect the quantum transition at n s k y /n e i ~ g even for disordered arrangement of skyrmions. 

In conslusion, we have investigated many-skyrmion ground state at low skyrmion concentration. We have shown 
that the skyrmions are bound in pairs, those forming triangular Wigner lattice. The quantun phase transition occurs 
in this lattice at certain skyrmion concentration. 
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FIG. 1. Binding of two skyrmions. Two skyrmion energy as a function of their separation has a minimum at zero separation. 

FIG. 2. Two phases of skyrmion lattice. In "insulating" phase, skyrmion phase strongly fluctuates so that x, y spin components 
have zero average. In "superconducting" phase, skyrmion phases are fixed and antiferromagnetically arranged 

FIG. 3. The quantum phase transition. Critical skyrmion density is plotted in the region of realistic g. Above the line, the 
system is in "superconducting" state. Numbers denote the equilibrium spin of isolated skyrmion pairs at given g. The crtical 
line drops to zero any time the spin state changes. 



4 



